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Abstract 



We derive a mesoscopic description of the behavior of a simple financial market where 
the agents can create their own portfolio between two investment alternatives: a stock 
and a bond. The model is derived starting from the Levy-Levy- Solomon microscopic 
model [141 115j using the methods of kinetic theory and consists of a linear Boltzmann 
equation for the wealth distribution of the agents coupled with an equation for the 
price of the stock. From this model, under a suitable scaling, we derive a Fokker-Planck 
equation and show that the equation admits a self-similar lognormal behavior. Several 
. numerical examples are also reported to validate our analysis. 



Keywords: wealth distribution, power-law tails, stock market, self-similarity, kinetic equa- 
tions. 
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-sj- ! 1 Introduction 

t> 

In recent years, physicists have been growing more and more interested in new interdis- 
ciplinary areas such as sociology and economics, originating what is today named socio- 
economical physics [H [3j 01 El [13l Q31 [201 12H ES] • This new area in physics borrows several 
methods and tools from classical statistical mechanics, where complex behavior arises from 
relatively simple rules due to the interaction of a large number of components. The moti- 
vation behind this is the attempt to identify and characterize universal and non-universal 
features in economical data in general. 

A large part of the research in this area is concerned with power-law tails with universal 
exponents, as was predicted more than one century ago by Pareto [H I16| [23] . In particular, 
by identifying the wealth in an economic system with the energy of a physical system, the 
application of statistical physics makes it possible to understand better the development of 
tails in wealth distributions. Starting from the microscopic dynamics, mesoscopic models 
can be derived with the tools of classical kinetic theory of fluids [H El [131 ED [22J [25] . 
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In contrast with microscopic dynamics, where behavior often can be studied only em- 
pirically through computer simulations, kinetic models based on PDEs allow us to derive 
analytically general information on the model and its asymptotic behavior. For example, 
the knowledge of the large-wealth behavior is of primary importance, since it determines a 
posteriori whether the model can fit data of real economies. 

In some recent papers, the explicit emergence of power laws in the wealth distribution, 
with Pareto index strictly larger than one, has been proved for open market economies where 
agents can interact through binary exchanges together with a simple source of speculative 
trading [U E2 E5] • 

The present work is motivated by the necessity to have a more realistic description of the 
speculative dynamics in the above models. To this end, we derive a mesoscopic description 
of the behavior of a simple financial market where a population of homogeneous agents 
can create their own portfolio between two investment alternatives: a stock and a bond. 
The model is closely related to the Levy-Levy-Solomon (LLS) microscopic model in finance 
[14} 1 15 j . This model attempted to construct from simple rules complex behavior that could 
then mimic the market and explain the price formation mechanism. As a first step towards 
a more realistic description, we derive and analyze the model in the case of a single stock 
and under the assumption that the optimal proportion of investments is a function of the 
price only. In principle, several generalizations are possible (different stocks, heterogeneous 
agents, a time-dependent optimal proportion of investments, . . .), and we leave them for 
future investigations. 

In our non-stationary financial market model, the average wealth is not conserved and 
this produces price variations. Let us point out that, even if the model is linear since no 
binary interaction dynamic between agents is present, the study of the large time behavior 
is not immediate. In fact, despite conservation of the total number of agents, we don't have 
any other additional conservation equation or entropy dissipation. Although we prove that 
the moments do not grow more than exponentially, the determination of an explicit form 
of the asymptotic wealth distribution of the kinetic equation remains difficult and requires 
the use of suitable numerical methods. 

A complementary method to extract information on the tails is linked to the possibility 
to obtain particular asymptotics which maintain the characteristics of the solution to the 
original problem for large times. Following the analysis in [5J, we shall prove that the 
Boltzmann model converges in a suitable asymptotic limit towards a convection-diffusion 
equation of Fokker-Planck type for the distribution of wealth among individuals. Other 
Fokker-Planck equations were obtained using different approaches in [H [2TJ, [18] . 

In this case, however, we can show that the Fokker-Planck equation admits self-similar 
solutions that can be computed explicitly and which are lognormal distributions. One is then 
led to the conclusion that the formation of Pareto tails in the wealth distribution observed 
in [U [5] is a consequence of the interplay between the conservative binary exchanges having 
the effect of redistributing wealth among agents and the speculative trading causing the 
growth of mean wealth and social inequalities. 

The rest of the paper is organized as follows. In the next section, we introduce briefly 
the microscopic dynamic of the LLS model. The mesoscopic model is then derived in 
Section 3 and its properties discussed in Section 4. These properties justify the asymptotic 
procedures performed in section 5. The model behavior together with its asymptotic limit 
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is illustrated by several numerical results in section 6. Some conclusions and remarks on 
future developments are then made in the last section. 



2 The microscopic dynamic 

Let us consider a set of financial agents % = 1, . . . , N who can create their own portfolio 
between two alternative investments: a stock and a bond. We denote by Wi the wealth of 
agent i and by raj the number of stocks of the agent. Additionally we use the notations S 
for the price of the stock and n for the total number of stocks. 

The essence of the dynamic is the choice of the agent's portfolio. More precisely, at each 
time step each agent selects which fraction of wealth to invest in bonds and which fraction 
in stocks. We indicate with r the (constant) interest rate of bonds. The bond is assumed 
to be a risk-less asset yielding a return at the end of each time period. The stock is a risky 
asset with overall returns rate x composed of two elements: a capital gain or loss and the 
distribution of dividends. 

To simplify the notation, let us neglect for the moment the effects due to the stochastic 
nature of the process, the presence of dividends, and so on. Thus, if an agent has invested 
"fiWi of its wealth in stocks and (1 — ^i)wi of its wealth in bonds, at the next time step in 
the dynamic he will achieve the new wealth value 

^ = (1-7^(1 + r)+ 7i ti; i (l + x), (1) 

where the rate of return of the stock is given by 

S' -S 



S ' 



(2) 



and S' is the new price of the stock. 
Since we have the identity 

we can also write 



'jiWi = mS, (3) 



w\ = Wi + Wi(l - %)r + wai ( 4 ) 
= Wi + (wi - UiS)r + m(S' - S). (5) 

Note that, independently of the number of stocks of the agent at the next time level, it is 
only the price variation of the stock (which is unknown) that characterizes the gain or loss 
of the agent on the stock market at this stage. 

The dynamic now is based on the agent choice of the new fraction of wealth he wants 
to invest in stocks at the next stage. Each investor i is confronted with a decision where 
the outcome is uncertain: which is the new optimal fraction j'- of wealth to invest in stock? 
According to the standard theory of investment each investor is characterized by a utility 
function (of its wealth) U(w) that reflects the personal risk taking preference |12| . The 
optimal 7^ is the one that maximizes the expected value of U(w). 
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Different models can be used for this (see |X5|, 129]). for example, maximizing a von 
Neumann-Morgenstern utility function with a constant risk aversion of the type 

U(w) = (6) 
1 — a 

where a is the risk aversion parameter, or a logarithmic utility function 

U{w) = log(w). (7) 

As they don't know the future stock price S", the investors estimate the stock's next 
period return distribution and find an optimal mix of the stock and the bond that maximizes 
their expected utility E[U]. In practice, for any hypothetical price S h , each investor finds the 
hypothetical optimal proportion ^{S h ) which maximizes his/her expected utility evaluated 
at 

v$(S h ) = (1 - ^)w'(l + r) + tM(1 + AS h )), (8) 

where x'(S h ) = (S h — S')/S' and S' is estimated in some way. For example in [15] the 
investors expectations for x' are based on extrapolating the past values. 

Note that, if we assume that all investors have the same risk aversion a in ([6]), then 
they will have the same proportion of investment in stocks regardless of their wealth, thus 
1 h { S h )= 1 h (S h )- 

Once each investor decides on the hypothetical optimal proportion of wealth 7^ that 
he/she wishes to invest in stocks, one can derive the number of stocks n^(S h ) he/she wishes 
to hold corresponding to each hypothetical stock price S h . Since the total number of shares 
in the market n, is fixed there is a particular value of the price S' for which the sum of 
the n^(S h ) equals n. This value S' is the new market equilibrium price and the optimal 
proportion of wealth is 7^ = j^(S'). 

More precisely, following [15], each agent formulates a demand curve 



characterizing the desired number of stocks as a function of the hypothetical stock price S h . 
This number of share demands is a monotonically decreasing function of the hypothetical 
price S h . As the total number of stocks 

JV 

n 

i=i 



E n * ( 9 ) 



is preserved, the new price of the stock at the next time level is given by the so-called 
market clearance condition. Thus the new stock price S' is the unique price at which the 
total demand equals the supply 

N 

Y^nKS')=n. (10) 

i=l 

This will fix the value w' in (pQ) and the model can be advanced to the next time level. To 
make the model more realistic, typically a source of stochastic noise, which characterizes all 
factors causing the investor to deviate from his/her optimal portfolio, is introduced in the 
proportion of investments 7« and in the rate of return of the stock x 1 . 
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3 Kinetic modelling 



We define / = f(w,t), w € M + , t > the distribution of wealth w, which represents the 
probability for an agent to have a wealth w. We assume that at time t the percentage of 
wealth invested is of the form 7(£) = /J,(S) +£, where £ is a random variable in [—z, z], and 
z = min{— n(S), 1 — fJ>(S)} is distributed according to some probability density 3>(//(,f>), £) 
with zero mean and variance Q 2 . This probability density characterizes the individual 
strategy of an agent around the optimal choice (i(S). We assume $ to be independent of 
the wealth of the agent. Here, the optimal demand curve //(•) is assumed to be a given 
monotonically non- increasing function of the price S > such that < /i(0) < 1. 

Note that given f(w,t) the actual stock price S satisfies the demand-supply relation 

S=-E[ 7 w], (11) 
n 

where E[X] denotes the mathematical expectation of the random variable X and f(w,t) 
has been normalized 

poo 

f(w, t)dw = 1. 





More precisely, since 7 and w are independent, at each time t, the price S(t) satisfies (see 
Figure [I} 

S(t) = -E[j\E[w] = -^(S(t))w(t), (12) 
n n 

with 



w(t):=E[w]= f(w,t)wdw (13) 
Jo 

being the mean wealth and by construction, 

At the next round in the market, the new wealth of the investor will depend on the 
future price S' and the percentage 7 of wealth invested according to 

w'(S', 7, t?) = (1 - 7)w(l + r) + jw(l + x{S', 77)), (14) 

where the expected rate of return of stocks is given by 

In the above relation, D > represents a constant dividend paid by the company and 77 
is a random variable distributed according to Q(rj) with zero mean and variance a 2 , which 
takes into account fluctuations due to price uncertainty and dividends [151 E] • We assume 
r] to take values in [— d, d] with < d < S' + D so that w' > and thus negative wealths are 
not allowed in the model. Note that equation (|15p requires estimation of the future price 
S', which is unknown. 
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Figure 1: Example of equilibrium price 



The dynamic is then determined by the agent's new fraction of wealth invested in stocks, 
7'(£') = l^(S') + where £' is a random variable in [—z 1 , z'\ and z' = min{/i(5'), 1 — n(S')} 
is distributed according to <E>(/i(S"), £'). We have the demand-supply relation 



n 

which permits us to write the following equation for the future price: 

S' = ^E[j']E[w'] = \i{S')E[w'}. 



Now 
thus 



w'(S', 7, rj) = w(l + r) + r yw(x(S' , rj) — r), 



E[w'} = E[w](l+r) + E[jw](E[x(S',ri)]-r) 

'S'-S + D 



rE(t)(l + r) + n(S)w(t) 



S 



This gives the identity 



1 



S' = -fi(S')w{t) 
n 



(l + r) + M (5) 



g - S + £> 
5 



(16) 

(17) 
(18) 

(19) 
(20) 

(21) 



Using equation (|12p we can eliminate the dependence on the mean wealth and write 

M (S') 



5' 



[<il-^S))S(l + r) + ^S)(S' + D)] 



(i-»(s)Ms>) 
(i-^s')MS) { + ' + i 



M (5') 



M (5')' 



(22) 
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Remark 3.1 

The equation for the future price deserves some remarks. 



Equation ([22]) determines implicitly the future value of the stock price. Let us set 

9{s) = ^r 5 - 

Then the future price is given by the equation 

g(S')=g(S)(l + r) + D 

for a given S. Note that 

dg(S) _ _MS)_J^_ 1 - mQS) 

so the function 5(5) is strictly increasing with respect to S. This guarantees the 
existence of a unique solution 

S' = g- 1 (g(S)(l+r)+D)> S. (23) 

Moreover, if r = and D = 0, the unique solution is S' = S and the price remains 
unchanged in time. 

For the average stock return, we have 

1 j (i-h(s')ms) + s(i-^(S')y [M) 

where 

x(S') = E[x(S', V )] = S '~ S s + D . (25) 

Now the right hand side of (f2H| has non-constant sign since fJ>(S') < (i(S). In partic- 
ular, the average stock return is above the bonds rate r only if the (negative) rate of 
variation of the investments is above a certain threshold 

»(S>)-ti(S)_ s> D 



M(S>(S') " (l + r) 

In the constant investment case yu(-) = C, with C € (0, 1) constant, then we have 
g (S) = (1 - C)S/C and 

S' = (l + r)S + T ^D, 

which corresponds to a dynamic of growth of the prices at rate r. As a consequence, 
the average stock return is always larger then the constant return of bonds: 
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By standard methods of kinetic theory [2] , the microscopic dynamics of agents originate 
the following linear kinetic equation for the evolution of the wealth distribution 



df(w,t) r d r* / ., \ i 



dt 



(p(!w -> w)—^f('w,t) - P(w «/)/(«;, t) J de^. (26) 



The above equation takes into account all possible variations that can occur to the distri- 
bution of a given wealth w. The first part of the integral on the right hand side takes into 
account all possible gains of the test wealth w coming from a pre-trading wealth 'w. The 
function {3('w — > w) gives the probability per unit time of this process. 
Thus 'w is obtained simply by inverting the dynamics to get 



'w 



w 



, j(S,r,,t) = l + r + 7 (£) (*(£?', t?) -r), (27) 



where the value 5" is given as the unique fixed point of (|17p . 

The presence of the term j in the integral is needed in order to preserve the total number 
of agents 

dtl /K*) d ™ = - 

The second part of the integral on the right hand side of (|26p is a negative term that takes 
into account all possible losses of wealth it; as a consequence of the direct dynamic (|14p , the 
rate of this process now being f3(w — > w'). In our case, the kernel f3 takes the form 

P(w^w') = $(p(S),£)e(ri). (28) 

The distribution function <&(//(S), £), together with the function //(■), characterizes the 
behavior of the agents on the market (more precisely, they characterize the way the agents 
invest their wealth as a function of the actual price of the stock). 

Remark 3.2 In the derivation of the kinetic equation, we assumed for simplicity that the 
actual demand curve //(•) which gives the optimal proportion of investments is a function 
of the price only. In reality, the demand curve should change at each market iteration 
and should thus depend also on time. In the general case where each agent has a wealth- 
dependent individual strategy, one should consider the distribution f(j, w, t) of agents hav- 
ing a fraction 7 of their wealth w invested in stocks. 



4 Properties of the kinetic equation 

We will start our analysis by introducing some notations. Let A^o be the space of all 
probability measures on R + and by 

M p = \^>eM :j |i9| p ^(t?)^ < +oo,p > o| , (29) 

we mean the space of all Borel probability measures of finite momentum of order p, equipped 
with the topology of the weak convergence of the measures. 
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Let J 7 P (]R + ), p > 1 be the class of all real functions on M + such that g(0) = g'(0) = 
and </ m ) (v) is Holder continuous of order 5, 

ii (m)ll |g (m) W-g (m) HI , ,„ n s 
\\g ( ' l> \\s = sup : -= < co, (30) 

v^w \V-W\ d 

the integer m and the number < 5 < 1 be such that m + 5 = p, and </ m ) denote the m-th 
derivative of g. 

Clearly the symmetric probability density which characterizes the stock returns be- 
longs to M. v for all p > since 

rd 

\r)\ p Q(r))dri < \d\ p . 



Moreover, to simplify computations, we assume that this density is obtained from a given 
random variable Y with zero mean and unit variance. Thus of variance a 2 is the density 
of crY. By this assumption, we can easily obtain the dependence on a of the moments of 
0. In fact, for any p > 2, 



f \rj\ p @(rj)dri = E (\aY\ p ) = a p E (\Y\ P ) 
J-d 

Note that equation (|26[) in weak form takes the simpler form 



j poo foo fD 



dt 



-D J- 



- f{w,t)^{w)dw= / / ^( f M(S),0&(r ] )f(w,t)(H^)-Hw))^dr ] dw. (31) 



By a weak solution of the initial value problem for equation (|26p corresponding to the 
initial probability density fo(w) G Ai p , p > 1, we shall mean any probability density 
/ E C 1 (M+, satisfying the weak form (|3ip for t > and all <fi € ^(K^.), and such that 
for all ^J P (R + ), 

/■oo rco 

lim / f(w,t)(j)(w)dw= / fo(w)(j}(w) dw. (32) 
Jo Jo 

The form (I31D is easier to handle, and it is the starting point to study the evolution of 
macroscopic quantities (moments). The existence of a weak solution to equation (j26|) can 
be seen easily using the same methods available for the linear Boltzmann equation (see |28j 
and the references therein for example). 

From (13ip follows the conservation of the total number of investors if (f>(w) = 1. The 
choice <j)(w) = w is of particular interest since it gives the time evolution of the average 
wealth which characterizes the price behavior. In fact, the mean wealth is not conserved in 
the model since we have 

d f°° ( ( S' — S + D \\ f°° 

^ J f(w,t)wdw=\r + fi(S){ r\ \ J f(w,t)wdw. (33) 

Note that since the sign of the right hand side is nonnegative, the mean wealth is nonde- 
creasing in time. In particular, we can rewrite the equation as 

jw(t) = ((1 - p(S))r + KS)x(S')) w(t). (34) 



9 



From this we get the equation for the price 
where S' is given by (I22p and 

Now since from (I24p it follows by the monotonicity of that 

x(S') < M d Mr + D 



S(0)(l-/z(S(0)))' 
using (|3~4"|) we have the bound 

w(f) < u)(0) exp (Mi) . (36) 
From (|12p we obtain immediately 

< .V \. exp (Mt) , 



M (S(t)) " /x(5(0)) 
which gives 

5(t) < 5(0) exp (Mi) . (37) 



Remark 4.1 For a constant fi(-) = C, C £ (0, 1) we have the explicit expression for the 
growth of the wealth (and consequently of the price) 

nD 

w(t) = w(0) exp(ri) - (1 - exp(rt)) — . (38) 

1 — 



Analogous bounds to (J36J) for moments of higher order can be obtained in a similar 
way. Let us consider the case of moments of order p > 2, which we will need in the sequel. 
Taking <p{w) = w p , we get 

j poo poo pd pz 

- wPf(w,t)dw= / / ^(S),OQ(ri)f(w,t)(w' p -wndidr,dw. (39) 

J " J — d J — z 

Moreover, we can write 

w'p = w p +pw p ~ 1 (w' -w) + ^p(p - l)w p ' 2 (w' - w) 2 , 
where, for some < $ < 1, 

w = $w' + (1 — $)w. 

Hence, 

roo pd pz 

/ / <S>(n(S),£)G(ri)f(w,t)(w' p - w p )d£dr]dw 

/0 J-dJ-z 
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roc rd 
'0 J-d 



From 



$(n(S),£)@{v)f(w,t)(pw p - 1 (w' -w) + -p{p-l)w p - 2 (w' - wf)d£ dv dw 

r°° i 
= p((l - /i(5))r + / i^/(ti>, *) + - 1) 

■/ o z 

oo /»c£ /-Z 

/ / ^(fi(S),^)e(r])f(w,t)w p - 2 w 2 ( y (l -j)r + jx(S',n)) 2 d£dr]dw. 

J-dJ-z 



w p - 2 = w p - 2 (l + - 7 )r + 7x(5 / , 7?))) p ~ 2 < w p ~ 2 (l + r + |x(S', 7?)|) p - 2 



< C p w; p - 2 (1 + r p ~ 2 + \x(S',r))\ p - 2 ) 

and 

((1 _ 7 ) r + 1X ( S >, v ))) 2 < 2(r 2 + x(S', r?) 2 ), 

we have 

"oo rd rz 

$(n(S),€)G(rj)f(w, t)w p - 2 w 2 {{\ - j)r + jx(S', rj)) 2 d£ dq dw 



Since 



d J — z 
•oo rd 

<2C P I / Q(r])f(w,t)w p (l + r p - 2 + \x(S',r])\ p - 2 )(r 2 + x(S',r]) 2 )dr]dw. 
J-d 



J_ Q(r,)\x(S , ,r,)Pdr, < -| ((5' - Sf + D p + <j p E{\Y\ p )) , (40) 
we finally obtain the bound 



-/ w p f(w,t)dw < A P (S) I w p f(w,t)dw, (41) 



where 

A P (S) = p((l-fi(S))r + ^S)x(S')) 

+ Pip ~ l)C p [r p + (1 + r p - 2 ) (l + - S) 2 + D 2 + ( ) N ) 

+ ^((S'-S) p + D p + a p E(\Y\ p )) 

and C p , c p , c p _2 and C2 are suitable constants. 
We can summarize our results in the following 

Theorem 4.1 Let the probability density /o S Mp, where p = 2 + 5 for some 8 > 0. Then 
the average wealth is increasing exponentially with time following A3 6)). As a consequence, 
if n is a non- increasing function of S, the price does not grow more than exponentially as 
in |ff7| ). Similarly, higher order moments do not increase more than exponentially, and we 
have the bound 
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5 Fokker-Planck asymptotics and self-similar solution 



The previous analysis shows that in general it is difficult to study in detail the asymptotic 
behavior of the system. In addition, we must take into account the exponential growth of the 
average wealth. In this case, one way to get information on the properties of the solution 
for large time relies on a suitable scaling of the solution. As is usual in kinetic theory, 
however, particular asymptotics of the equation result in simplified models (generally of 
Fokker-Planck type) whose behavior is easier to analyze. Here, following the analysis in 
[5j [22] and inspired by similar asymptotic limits for inelastic gases [H [25] , we consider the 
limit of large times in which the market originates a very small exchange of wealth (small 
rates of return r and x). 

In order to study the asymptotic behavior of the distribution function f(w,t), we start 
from the weak form of the kinetic equation 



r'OO pd 

f(w,t)cf>(w)dw = I / / ^(S)^)^)!^,^^') - <f>(w))d^dr,dw (42) 



d_ 

dt Jo JO J-dJ - 

and consider a second-order Taylor expansion of 4> around w, 

4>{w') — 4>{w) = w(r + j(x(S', rj) — r))(f>'(w) + -w 2 (r + r y(x(S' , rf) — r)) 2 cp" (w) , 
where, for some < $ < 1, 

w = fin)' + (1 — €)w. 
Inserting this expansion into the collision operator, we get 

rd 



z 

d 

$(/u(5), £)e(rj)f(w, t)w{r + j(x(S', rj) - r))<f>'(w)d£ drj dw 



1 



2 



oo I'd 



-dJ- 



+- / / / ^(S),()e(r])f(w,t)w 2 (r + 7 (x(5',77) - r)) 2 <p" (w)d£ dr] dw 



+R r (S, S'), 



where 



-i roo rd pz 

Rr(S,S') = - *MS),Oe(r,)f(w,t) 

z JO J-dJ-z ^g-j 

w 2 (r + j(x(S', rj) - r)) 2 {<j)"{w) - <p"(w)) d£ dr) dw. 

Recalling that E[£] = 0, E[rj\ = 0, E[£ 2 ] = C 2 and E[rf] = a 2 , we can simplify the above 
expression to obtain 

oo rd fz 

$(v(S),0<3>(r])f(w,tmw') - <t>(w))dtdr]d W 



J-d 
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f(w, t)w + fi(S) ^ — r^j^j (j/(w) dw 

1 f°° , f \ 1 { 2 sn^U(S'-S) 2 a 2 + D 2 „^S f -S 

+ - / /(«;, t)w 2 r 2 + (C 2 + MS) 2 ) V a9 ^ + a9 + 2D- 



2Jo V V S 2 S 2 S 2 

2 S'-S + D\ fff-S + D \\ ... , 
+ r 2 - 2r J + 2r / u(S) f r ^'(u>) dio 

+i? r (,S, 5'). 
Now we set 

r = ri, f(w,r)=f(w,t), S(r) = S(t), fi(S) = fi(S), 

which implies that f(w, r) satisfies the equation 

d f°° ~ 

J f(w,T)(j)(w)dw 

' f(w, t)w [l + /i(5) + ^~ S -lj^j 4>'{w) dw 
+- / f(w, t)w 2 [r + (C 2 + /i(S) 2 ) ^ + ^— + 2L>- 



2 Jo V V ^<S* 2 rS 2 rS 2 

+-^(5,5'). 

r 

Now we consider the limit of very small values of the constant rate r. In order for such a 
limit to make sense and preserve the characteristics of the model, we must assume that 

lim — = v, lim — = A. (44) 
r->0 r r— >o r 

First let us note that the above limits in (j22[) imply immediately that 

lim 5' = S. (45) 

r— >0 

We begin by showing that the remainder is small for small values of r. 
Since <p £ J r 2+ ( 5(M+) and [u> — u)j = i9\w' — w\, 

\4>"{w) - <j>"{w)\ < U"\\s\w - w\ 5 < U"\\s\w' - w\ 5 . (46) 

Hence 

\-R r (s,s')\ < ^ff f^WUjeW 

r 2r J J_ d J_ z 

|(1 - j)r + 7 x(5', r?) | 2+5 /(^, r)w 2+5 (i£ cfy dw. 

By the inequality 

|(1 _ 7 ) r + 7 x(5',r ? )| 2+5 < 2 2 + 5 (r 2+5 + |x(S',7/)| 2+ *) 
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and (|4"0j) . we get 



\-R r (S,S')\ < 2 1+s U" \\ s - 
r 



r l+5 + C2+6 



((§' - S) 2+s + D 2+s + <j 2+s E(\Y\ 2+6 )) ) / f(w,T)w 2+6 dw. 



„g2+5 



2+S\ 



.,2+5, 



As a consequence of (}4"4"]) - ([4"5j) . from this inequality it follows that R r (S,S') converges to 
zero as r — > if 

/•oo 

,,2+5 



W 2+d f( W ,T)dw 

remains bounded at any fixed time r > 0, provided the same bound holds at time r = 0. 
This is guaranteed by inequality (jlTj) since A P (S) — >■ in the asymptotic limit defined by 
(31. 

Next we write 

£(£') = + (5' - S)fl(S) + 0((S' - S) 2 ), 

where 

fcS) = «if> < 0. 

(ZD 



(47) 



Then, using the above expansion from (1441) in (j22H . we obtain 



lim 

r— >o r 



k(S) 5 + 



^A 



1-£(S) 



with 



< k{S): 



def fi(S)(l - fl(S)) 



< 1. 



£(S)(1- £(5)) 

Now, sending r — >• under the same assumptions, we obtain the weak form 

d ro ° 



d,T 



f(w, r)(j)(w)dw 



(48) 



£(S)(k(S) - 1) + 1 A 



i - 



s 



/(w, r)w<j) (w) dw 



HK§) 2 + C 2 



s „, ^ y f(w,r)w 4> {w)dw. 
This corresponds to the weak form of the Fokker-Planck equation 



or equivalently 



— / — 

dr dw 



-A(r)wf + ^B(r)^-w 2 f 



(49) 
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with 



Mr) = l + A(g)fwS)-l) + A(S)( ;^ ) 1) + 1 |) (50) 

m _ (i(sr±n v . " (S1) 

Thus we have proved 

Theorem 5.1 Lei t/ie probability density /q € -M p; where p = 2 + 5 for some 5 > 0. Then, 
as r — > 0, a — > 0, and D — > in such a way that a 2 = vr and D = Xr, the weak solution to 
the Boltzmann equation l3l\) for the scaled density f r (w,r) = f(v,t) with r = rt converges, 
up to extraction of a subsequence, to a probability density f(w,r). This density is a weak 
solution of the Fokker-Planck equation 

We remark that even for the Fokker-Planck model, the mean wealth is increasing with 
time. A simple computation shows that 

d /*°° /'oc 
w(t) = — / f (w , t)w dw = A(t) / f(w, t)w dw = A{t)w{t). (52) 
dr Jo Jo 

Using dM|), we get the bound 

(1 - KS))w(t) + n\< w(r) < w{t) + (53) 

1 - y.{S) 

Similarly, for the second-order moment we have 

.1 roc poo 

g(r) = — / f(w, t)w 2 dw = (2A(t) + B(t)) / f(w, t)w 2 dw = (2A(t) + B(r))e(r). (54) 
dr Jo Jo 

In order to search for self-similar solutions, we consider the scaling 

f(w,r) = —g(x,T), x = log(w). 
w 

Simple computations show that g(x, T ) satisfies the linear convection-diffusion equation 
d , (B(t) \ d B(r) d 2 

which admits the self-similar solution (see [T7] for example) 

5fy r) ~ 1 exp ( (x + Kr)/2-a(r)) 2 \ 

9{X,T) ~ (2Kr)vr)V2 exp V 2 &(r) J' (55) 



where 

a(r) 



f'T f'T 

I A(s)ds + C 1 , b(r)= B(s)ds + C 2 . 
Jo Jo 



15 



Reverting to the original variables, we obtain the lognormal asymptotic behavior of the 
model, 

The constants C\ = a(0) and Ci = 6(0) can be determined from the initial data at t = 0. 
If we denote by w{0) and e(0) the initial values of the first two central moments, we get 

e(0) 



d = log(«J(0)), C 2 = log 



(tD(O)) 



Finally, a direct computation shows that 



a(r) = f ^ ds + C i = logHT)) (57) 
Jo w{s) 

and 

Jo w(s)J \{w{t)) 2 J 

Remark 5.1 

• If we assume £ and a are of the same order of magnitude, in the Fokker-Planck limit 
the noise introduced by the agents' deviations with respect to their optimal behavior 
does not play any role and the only source of diffusion is due to the stochastic nature 
of the returns. 

• In the case of constant investments ) = C, C € (0, 1) we have the simplified 
Fokker-Planck equation 



— / = — 

dr dw 



C A\ ~ 1(C 2 + C 2 ) d 2 ~ 
1 H — — wf + _ v—w l $ 



l-CS 2 S 2 dw 



It is easy to verify that for such a simple situation, the pair of ordinary differential 
equations for the evolution of the first two central moments, (|52[) and (|54|) . can be 
solved explicitly. 

Imposing the conservation of the mean wealth with the scaling 

/(-.r) = ^/(.,r), (59) 

W(T) W{T) 

we have the diffusion equation 

3 t . , B(t) d 2 2 1 



This yields the asymptotic lognormal behavior 



f (v T) = I exp f Q0Z{v)+l0Z(VW)/m)) 2 ) m 

JK,> i;(21og( J B(r)/^(0) 2 )7r) 1 /2 p ^ 2 log(E( T )/w(Q) 2 ) J 

with 

f(v,r)vdv = w(0), E(t) = \ f(v,r)v 2 dv. 
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6 Numerical examples 



In this section we report the results of different numerical simulations of the proposed 
kinetic equations. In all the numerical tests, we use N = 1000 agents and n = 10000 shares. 
Initially, each investor has a total wealth of 1000 composed of 10 shares, at a value of 50 per 
share, and 500 in bonds. The random variables £ and r] are assumed distributed according 
to truncated normal distributions so that negative wealth values are avoided (no borrowing 
and no short selling). In Tests 1 and 2 we compare the results obtained with the Monte 
Carlo simulation of the kinetic model to a direct solution of the price equation f)35[) . In the 
last test case we consider the time-averaged Monte Carlo asymptotic behavior of the kinetic 
model and compare its numerical self-similar solution with the explicit one computed in the 
last section using the Fokker-Planck model. 

Test 1 

In the first test we take a riskless interest rate r = 0.01 and an average dividend growth rate 
D = 0.015 and assume that the agents simply follow a constant investments rule, /i(-) = C, 
with C € (0, 1) constant. As a consequence of our choice of parameters we have C = 0.5 and 
the evolution of the mean wealth and of the price in the kinetic model are known explicitly 
(|38p . We report the results after 400 stock market iterations with £ and rj/S(0) distributed 
with standard deviation 0.2 and 0.3 respectively. In Figure [2] we report the simulated price 
behavior together with the evolution computed from (|35p . The fraction of investments in 
time during the Monte Carlo simulation fluctuates around its optimal value and is given in 
Figure El 

Test 2 

In the next test case we take the same parameters as in Test 1 but with a non constant 
profile /i(-). More precisely we take a monotone decreasing exponential law 

li(S) = Ci + (1 - d)e- C2S 

with C\ = 0.2 and Ci = log((l — Ci)/(0.5 — C\))/So m 0.02 so that the price equation is 
satisfied for Sq = 50. We have 0.2 < //(•) < 0.5. The results for the price evolution and 
the investments behavior are plotted in Figures 0] and The solution for the price in the 
kinetic equation has been computed by direct numerical discretization of (|35p . Note that 
the final price is approximatively 5 times smaller then the one in the constant investment 
case with [x = 0.5. In Figure El we compare the behavior of the mean wealth in Tests 1 and 
2 to the exponential growth at a rate r obtained with simple investments in bonds. We can 
observe that the time decay of investments in Test 2 is fast enough to produce a wealth 
growth below the rate r. On the contrary, as observed in Section 3, a constant investment 
strategy produces a curve above this rate. Finally, in Figure [71 we plot the averaged wealth 
distributions at the final computation time on a log- log scale together with a lognormal fit. 
The results show lognormal behavior of the tails even for the Boltzmann model. Note that 
thanks to equation ([3]), the same distribution is observed for the number of stocks owned 
by the agents. 
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Figure 2: Test 1. Exponential growth of the price in time. Numerical simulation of the 
kinetic model. 
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Figure 3: Test 1. Fluctuations of the corresponding fraction of investments in time. 
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Figure 4: Test 2. Price behavior in time. Numerical simulation of the the kinetic model. 



0.25 



0.2 



0.15 




50 



Figure 5: Test 2. Fluctuations of the corresponding fraction of investments in time. 
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Figure 6: Tests 1 and 2. Behavior of the mean wealth in the kinetic model. The top curve 
refers to Test 1, the bottom curve to Test 2. The dashed line corresponds to the exponential 
growth at a rate equal to r. 




w 



Figure 7: Tests 1 and 2. Log-log plot of the mean wealth distribution together with a 
lognormal fitting. The top curve refers to Test 1, the bottom curve to Test 2. 
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Test 3 



In the last test case we consider the asymptotic limit of the Boltzmann model and compare 
its numerical self-similar solution with the explicit one computed in the last section using the 
Fokker-Planck model. To this end, we consider the self-similar scaling (|59p and compute the 
solution for the values r = 0.001, D = 0.0015 with £ and r]/S(0) distributed with standard 
deviation 0.05. We report the numerical solution for a constant value of /i = 0.5 at different 
times t = 50,200,500 in Figures [8] and A very good agreement between the Boltzmann 
and the lognormal Fokker-Planck solutions is observed, as expected from the results of the 
last section. We also compute the corresponding Lorentz curve L(F(w,t)) defined as 



The Gini coefficient is a measure of the inequality in the wealth distribution |10j . A value of 
corresponds to the line of perfect equality depicted in Figure [TUl together with the different 
Lorentz curves. It is clear that inequalities grow in time due to the speculative dynamics. 

7 Conclusions 

We have derived a simple linear mesoscopic model which describes a financial market under 
the assumption that the distribution of investments is known as a function of the price. The 
model is able to describe the exponential growth of the price of the stock and the growth 
of the wealth above the rate produced by simple investments in bonds. The long-time 
behavior of the model has been studied with the help of a Fokker-Planck approximation. 
The emergence of a power law tail for the wealth distribution of lognormal type has been 
proved. In order to produce the effect of a real financial market, with booms, cycles and 
crashes, the distribution of investments should be a function of time (the decision-making 
should be done by maximizing the expected utility) and one should consider heterogeneous 
populations of investors as in [14| 115] . In this case, the model should be modified and the 
time evolution of )J-(S, t) considered. Another interesting research direction is related to the 
possibility to introduce stock options into the model and to relate the kinetic approach to 
Black-Scholes type equations. All these subjects are actually under investigation and we 
hope to present other challenging results in the near future. 
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and the Gini coefficient G € [0, 1] 
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Figure 8: Test 3. Distribution function at t = 50, 200, 500. The continuous line is the 
lognormal Fokker-Planck solution. 
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w 

Figure 9: Test 3. Log-log plot of the distribution function at t = 50, 200, 500. The continu- 
ous line is the lognormal Fokker-Planck solution. 
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Figure 10: Test 3. The corresponding Lorentz curves. The Gini coefficients are G = 0.1, 
G = 0.2 and G = 0.3 respectively. 
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